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I.  INTRODUCTION 


1.1  Background 

The  need  to  develop  microelectronic  devices  capable  of  operating  at 
increased  performance  levels  with  high  reliability,  requires  better  methods  of 
thermal  control.  With  the  recent  trend  towards  increased  miniaturization  and 
component  density,  the  thermal  management  problem  increases  dramatically.  This 
is  particularly  true  in  the  microelectronic  components  proposed  for  use  in 
spacecraft  applications  where  the  heat  removal  capacity  of  the  thermal  control 
systems  governs  the  physical  size  of  both  the  individual  components  and  overall 
systems. 

In  1984,  Cotter  described  very  small,  micro  heat  pipes  for  cooling 
microelectronic  devices.  He  defined  a  micro  heat  pipe  as  "one  so  small  that  the 
mean  curvature  of  the  vapor-liquid  interface  is  necessarily  comparable  in 
magnitude  to  the  reciprocal  of  the  hydraulic  radius  of  the  total  flow  channel." 
Mathematically,  this  can  be  expressed  as 

1 

rc  « -  (1) 

rh 

where  is  the  mean  curvature  of  the  vapor-liquid  interface  and  r^  is  the 
hydraulic  radius  of  the  total  flow  channel.  This  relationship  can  be  better 
defined  by  nondimensionalizing  both  the  mean  curvature  of  the  vapor-liquid 
interface  and  and  the  hydraulic  radius,  i.e.,  dividing  both  by  rc,  prior  to 
substitution  into  equation  1.  This  yields  an  expression  which  more  precisely 
defines  a  micro  heat  pipe  as  one  which  satifies  the  condition  that 

—  >  1  (2) 

rh 
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In  practical  terms,  a  micro  heat  pipe  is  a  wickless,  noncircular  channel 
with  an  approximate  diameter  of  100  to  1000  ym  and  a  length  of  about  10  to  20 

mm.  Unlike  traditional  heat  pipes,  micro  heat  pipes  would  not  have  a  wick  but 

would  use  the  mean  curvature  of  the  channel  to  pump  the  liquid  from  the 

condenser  to  the  evaporator,  which  would  be  an  integral  part  of  the  device. 

Besides  the  difficulties  of  manufacturing,  cleaning,  and  filling  these  miniature 
heat  pipes,  they  would  be  very  sensitive  to  the  amount  of  working  fluid  present. 
Since  the  liquid-vapor  interface  changes  continually  along  the  pipe,  care  must 
be  taken  to  ensure  proper  wetted  conditions  without  flooding.  These  unique 
devices  may  play  an  important  role  in  future  microelectronic  devices  where 
thermal  loads  are  high  and  space  is  at  a  premium. 

1.2  Research  Objectives 

To  identify  and  better  understand  the  phenomena  which  govern  the 
performance  limitations  and  operating  characteristics  of  micro  heat  pipes,  a 
combined  experimental  and  analytical  investigation  was  conducted.  The 
objectives  of  this  investigation  were  twofold: 

i.  Develop  an  analytical  model  capable  of  predicting  the  steady-state 

operational  characteristics  and  performance  limitations  of  a  micro  heat 
pipe. 

ii.  Obtain  and  test  several  micro  heat  pipes  to  determine  the  effects  of 
the  amount  of  working  fluid,  the  shape  of  the  capillary  wicking 
structure,  the  operating  temperature,  the  evaporator  heat  flux,  and  the 
effects  of  gravity  on  the  heat  transport  capacity  and  performance 
limi tations . 

Presented  here  is  a  description  of  a  steady-state  analytical  model  and  the 
methodology  utilized  in  the  development  of  . that  model  along  with  a  description 
of  the  experimental  portion  of  the  investigation,  including  a  description  of  the 
test  article,  the  experimental  apparatus,  and  the  experimental  procedure 
followed.  In  Section  IV,  the  experimental  results  are  compared  with  the 
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predicted  results  obtained  from  the  steady-state  analytical  model. 

Figure  1  illustrates  the  general  shape  of  the  heat  pipe  being  modeled, 
along  with  the  liquid  configuration,  and  Fig.  2  illustrates  the  critical 
cross-sectional  dimensions.  As  shown,  the  outer  dimensions  are  1  mm  by  1  mm  and 
the  length  is  57  mm.  The  case  material  is  either  copper  or  silver  and  the 
working  fluid  is  water.  This  heat  pipe  is  considerably  smaller  than  most  heat 
pipes  and  because  of  this,  the  triangular  grooves  in  the  corners  of  the  heat 
pipe  can  serve  as  arteries  for  the  return  of  liquid  (i.e.,  the  wicking 
structure).  Because  of  the  size  of  the  heat  pipe  tested,  it  was  very  difficult 
to  measure  the  thermal-hydraulic  parameters  within  the  pipe  during  operation. 

For  this  reason,  in  addition  to  accomplishing  the  two  objectives  outlined  above, 
a  transient  analytical  model  was  developed  for  the  micro  heat  pipe.  This  model 
is  capable  of  predicting  the  transient  behavior  during  start-up  or  during 
variations  in  the  thermal  load  in  the  evaporator  portion  of  the  heat  pipe. 

Using  this  model,  the  time  required  to  reach  steady-state  can  be  predicted  along 
with  the  effect  of  variations  in  the  amount  and  type  of  fluid  and  the 
cross-sectional  shape.  A  detailed  description  of  the  methodology  used  in  the 
development  of  this  transient  model  along  with  a  comparison  of  the  steady-state 
experimental  results  and  the  steady-state  results  obtained  from  the  transient 
model  are  presented. 

II.  STEADY-STATE  MODEL 

Although  the  fundamental  phenomena  involved  in  the  operation  of  micro  heat 
pipes  may  be  different  from  those  occurring  in  larger  heat  pipes,  traditional 
steady-state  modeling  techniques  can  be  used  as  a  first  estimate  of  the 
operational  characteristics  and  performance  limitations.  Following  is  a 
discussion  of  the  traditional  limitations  of  heat  pipes  and  how  they  pertain  to 
the  configuration  described  above. 
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Fig.  1  Trapezodial  micro  heat  pipe. 
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Fig.  2  Cross-sectional  shape  and  dimensions  of  the  trapezoidal  micro  heat  pipe. 

2.1  The  Capillary  Limitation 

Because  of  the  size  and  wicking  structure  of  micro  heat  pipes,  we 
anticipated  that  they  will  be  limited  by  capillary  forces.  The  operation  of  a 
heat  pipe  requires  the  capillary  pumping  presure  to  be  greater  than  all  other 
resisting  pressures;  the  hydrostatic  pressure,  the  liquid  pressure,  and  the 
vapor  pressure.  The  governing  relation  is  expressed  by  the  equation 

GP  >  +AP  |  ,  +  AP,  +  AP  (3) 

c  -  |  |  1  v 

The  capillary  pumping  pressure  (APc)  is  the  pressure  difference  between  the 
evaporator  and  condenser  ends  of  the  heat  pipe  at  the  liquid  and  vapor 
interface.  The  hydrostatic  pressure  drop  (APii)  is  due  to  the  body  forces.  The 
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vapor  (AP^)  and  liquid  (AP^)  pressure  drops  are  due  to  viscous  forces  in  the 
vapor  and  liquid,  respectively. 

Two  methods  will  be  discussed,  both  of  which  describe  how  this  limit  can  be 
used  to  predict  steady-state  operation.  The  first  method  is  an  adapted  version 
of  conventional  heat  pipe  theory  while  the  second  method  has  been  specifically 
developed  for  micro  heat  pipes. 

2.1.1.  The  Conventional  Method.  The  conventional  method  involves  defining 
each  of  the  pressure  terms  in  equation  (1)  independently  and  then  combining  them 
to  solve  for  the  maximum  heat  flux. 

First,  it  is  necessary  to  develp  a  relationship  for  the  capillary  pumping 
pressure.  Reducing  the  LaPlace-Young  equation  results  in 


2  a 

2  a 

r 

ce 

r 

cc 

where  r  and  r  are  the  radii  of  curvature  of  the  fluid  at  the  evaporator  and 
condenser  regions,  respectively  and  cr  is  the  surface  tension  of  the  fluid. 
Because  of  vaporization  at  the  evaporator  region  of  the  heat  pipe,  the  liquid 
meniscus  recedes  in  the  triangular  groove.  Also,  in  the  condenser  section 
condensation  causes  the  liquid  meniscus  to  be  flooded.  This  resulting 
difference  in  radii  of  curvature  causes  a  pressure  difference  and  hence,  pumping 
of  the  liquid  from  the  condenser  to  the  evaporator.  It  is  generally  assumed 
that  rcc  approaches  infinity  during  steady-state  operation  which  allows  the 
capillary  pressure  to  be  only  a  function  of  the  evaporator  capillary  radius,  or 

AP  =  (2o/r  )  (5) 

c  ce  v 

The  wicking  structure  used  in  micro  heat  pipes  can  be  approximated  as  triangular 
grooves.  Chi  (1976)  found  that  for  triangular  grooves 
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=  v/cos(a) 


(6) 


r 

ce 

where  w  is  the  groove  width  and  a  is  the  half-included  angle. 

The  axial  hydrostatic  pressure  drop  is  a  result  of  the  body  forces 
acting  along  the  longitudinal  axis.  This  is  expressed  as 

AP  |  |  =  p1gLpsin(<f>)  (7) 

where  Lp  is  the  overall  length  of  the  heat  pipe  and  <|>  is  the  angle  the  heat 
pipe  makes  with  respect  to  the  horizontal. 

The  third  term  necessary  is  the  liquid  pressure  drop.  The  viscous 
forces  in  the  liquid  result  in  a  pressure  drop  which  resists  the  capillary 
flow  through  the  wick.  Since  the  liquid  pressure  gradient  may  vary  along 
the  longitudinal  axis  of  the  heat  pipe,  the  total  liquid  pressure  drop  can 
be  determined  by  integrating  the  pressure  gradient  over  the  length  of  the 
flow  passage,  or 


APX  (z) 


'z  dP. 

- dz 

0  dz 


(8) 


where  the  limits  of  integration  are  from  the  evaporator  end  to  the 
condenser  end  (z=0)  and  dP^/dz  is  the  gradient  of  the  liquid  pressure 
resulting  from  frictional  drag.  This  frictional  drag,  due  to  the  shear  stress, 
can  be  written  as 


dpx 

dz 


2t, 


(rh,l> 


(9) 


where  is  the  frictional  shear  stress  at  the  liquid-solid  interface  and 
r^  ^  is  the  hydraulic  radius,  defined  as  twice  the  cross-sectional  area 
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divided  by  the  wetted  perimeter. 

At  this  point,  the  Reynolds  number,  Re^,  and  drag  coefficient,  f^,  can  be 
introduced  as 


and 


Re!  = 


2(rh,l)plVl 


(10) 


f 


1 


2T-. 


plVl 


(ID 


where  is  the  local  liquid  velocity  which  is  related  to  the  local  heat 
flow  by  the  equation 


V1  = 


eViL 


(12) 


In  this  expression,  A  is  the  groove  cross-sectional  area,  e  is  the  wick 

S 

porosity,  and  L  is  the  latent  heat  of  vaporization. 

Combining  equations  (8)  through  (10)  with  (7)  yields, 


dp. 


dz 


{  (Vel)ul 


^eAg(rh,i)  Lpr 


This  equation  can  be  written  as 

dpx 

dz 


l  KALp  ) 
g  1 


where  K  is  the  permeability  expressed  as 


K  = 


(13) 


(14) 


(15) 


8 


For  laminar  flow,  (/^Re^)  is  constant  and  depends  only  on  the  passage  shape. 
Thus,  the  permeability,  K,  is  independent  of  the  flow  path.  For  triangular 
grooves 


K 


(16) 


where  d^  is  the  hydraulic  diameter  and  B  is  a  constant  with  a  value  of  30 
for  smooth  copper  or  silver. 

For  the  axial  direction  with  constant  heat  addition  and  removal, 
equation  (12)  can  be  substituted  into  equation  (6)  and  integrated  over  the 
length  of  the  heat  pipe  to  yield 


'  *1 

JCA  Lp. 
g  1 


(17) 


where  Le^.  is  the  effective  heat  pipe  length  defined  as 


L  =  0.5L  +  L  +  0.5L 

eff  e  a  c 


(18) 


Determination  of  the  vapor  pressure  drop  in  heat  pipes  is  complicated 
by  the  mass  addition  and  removal  in  the  evaporator  and  condenser 
respectively,  and  by  the  compressibility  of  the  vapor  phase.  Performing  a 
mass  balance  on  a  section  of  the  adiabatic  region  of  the  heat  pipe  ensures 
that  for  continued  operation,  the  liquid  mass  flow  rate  and  vapor  mass  flow  rate 
must  be  equal.  Because  of  the  large  difference  in  the  density  of 
these  two  phases,  the  velocity  of  the  vapor  phase  must  necessarily  be 
significantly  higher  than  the  velocity  of  the  liquid  phase.  For  this 
reason,  in  addition  to  the  pressure  gradient  resulting  from  frictional 
drag,  the  pressure  gradient  due  to  variations  in  the  dynamic  pressure  must 
also  be  considered.  Both  Chi  (1976)  and  Dunn  and  Reay  (1983)  have 
addressed  this  problem.  Chi  (1976)  found  that  upon  integration  of  the 
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vapor  pressure  gradient,  the  dynamic  pressure  effects  cancel.  The  result 
is  an  expression,  which  is  similar  to  that  developed  for  the  liquid, 


AP 

v 


'  C<*vRev>uv  1 


^<rh,v>  VvLj 


Le/A 


(19) 


where  r^  ^  is  the  hydraulic  radius  of  the  vapor  space  and  C  is  a  constant,  which 
depends  on  the  Mach  Number. 

As  mentioned  previously,  during  steady-state  operation,  at  any  axial 
position,  the  liquid  mass  flow  rate,  m^,  must  equal  the  vapor  mass  flow 
rate,  mv,  and  while  the  liquid  flow  regime  is  always  laminar,  the  vapor 
flow  may  be  either  laminar  or  turbulent.  It  is  therefore  necessary  to 
determine  the  vapor  flow  regime  as  a  function  of  the  heat  flux.  This  can 
be  accomplished  by  evaluating  the  local  axial  Reynolds  numher  in  the  vapor, 
defined  as 


Re 


v 


A  u  L 

v  \) 


(20) 


In  addition,  it  is  necessary  to  determine  if  the  flow  should  be  treated  as 
compressible  or  incompressible  by  evaluating  the  local  Mach  number,  defined  as 


Ma 


A  p  L(R  T  y) 
v  V  V  v 


.5 


(21) 


where  R^  is  the  gas  constant,  T^  is  the  vapor  temperature  and  y  is  the 
ratio  of  specific  heats,  which  is  equal  to  1.67,  1.4,  or  1.33  for 
monatomic,  diatomic,  and  polyatomic  vapor,  respectively  (Chi,  1976). 

Previous  investigations  summarized  by  Kraus  and  Bar-Cohen  (1983)  have 
demonstrated  that  the  following  combinations  of  these  conditions  can  be 
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used  with  reasonable  accuracy. 


Rev  <  2300,  Mav  <0.2 
(/vRev)  =  16 
C  =  1.00 

Re  <  2300,  Ma  <0.2 
v  -  v 

</vRev)  -  16 


Y  -  1' 

1  + 

Ma  2 

V 

2  > 

Re^  >  2300,  Ma^  <0.2 
(fv Rev)  =  0.038 


(22) 


(23) 


(24) 


Re^  >  2300,  Mav  >0.2 
(/v Rev)  =  0.038 


(25) 


Since  the  equations  used  to  evaluate  both  the  Reynolds  number  and  the 
Mach  number  are  functions  of  the  heat  transport  capacity,  it  is  necessary 
to  first  assume  the  conditions  of  the  vapor  flow.  Using  these  assumptions,  the 
maximum  heat  capacity,  Q  ,  can  be  determined  by  substituting  the 
values  of  the  individual  pressure  drops  into  equation  (1)  and  solving  for 
Qc  m-  The  overall  solution  for  the  heat  transport  capacity  as  limited  by 
the  capillary  forces  of  the  wicking  structure  is  then 
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-  -  PjgLsin(u) 
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ce 
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’  >-1  ' 

KL»1  + 


(26) 


with  the  values  for  C  and  /^Rev  found  iteratively,  depending  on  the  vapor 
flow  regime. 

2.1.2.  Cotter's  Method.  Cotter  (1984)  developed  a  steady-state  equation 
that  described  the  capillary  limit  specifically  for  micro  heat  pipes.  Since 
this  method  applies  only  to  micro  heat  pipes  many  simplifying  assumptions  were 
made.  These  assumptions  will  be  discussed  as  the  final  equation  is  developed. 

This  development  begins  with  restating  equation  (1).  By  summing  the 
pressures  throughout  the  fluid  flow  passage,  an  expression  for  the  liquid-vapor 
pressure  difference  can  be  obtained 


P 


v 


(27) 


where  and  P^  are  the  vapor  and  liquid  pressures  and  rm  is  the  mean  radius  of 
curvature  given  by 


1  1  1 
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r 

cc 


r 

ce 


(28) 


Simplifying  the  momentum  equation,  the  equation  of  the  fluid  in  steady-state 
parallel  flow  becomes 


U 


'6 

Sx 


+ 


dP 

dz 


(29) 


where  w  is  the  axial  fluid  velocity.  For  a  uniform  cross-sectional  area,  A, 
having  w=0  on  the  boundaries,  we  can  transform  to  new  dimensionless 
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variables  x'  ,  y'  ,  and  v'  given  by 


1/2  1/2  A  dP 

x  =  A  x' ;y  =  A  y'  ;v  = - w' 

y  dz 


Substituting  these  variables  into  equation  (29)  satisfies 


Sv '  Sw ' 

2~  +  2 

Sx'z  6y ' 2 


(30) 


(31) 


Again  substituting  the  dimensionless  variables  into  the  local  mass  flow 
equation, 


m 

=  P 


wdxdy 


the  result  is 

K'  PA2  dP 


8  ny  dz 


(32) 


(33) 


where  K'  is  a  shape  factor  given  by 


K'  =  8n 


w'dx'dy' 


(34) 


Values  for  K'  have  been  developed  by  Cotter  (1984)  and  are  presented  in 
Table  1. 


Table  1 

Flow  Shape  Factor  Values 


Shape  K' 


Circle 

1.000 

Regular  Hexagon 

0.964 

Square 

0.883 

Triangle,  60-60-60 

0.725 

Triangle,  45-55-90 

0.656 

Triangle,  30-60-90 

0.597 
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Solving  for  the  pressure  differential  in  both  the  vapor  and  liquid 
cases  and  substituting  this  into  equation  (27)  yields 


8nvJ”v  8nvm^  d  a 

K'  A2  K',  A2,  dz  R 

v  v  11 


(35) 


Realizing  that  for  steady-state,  the  local  liquid  and  vapor  mass  flow  rates 
are  equal  in  magnitude  and  proportional  to  the  local  heat  flux,  Q,  the  mass  flow 
rate  can  be  calculated  as 


m 


Q 


(z)  =  —  h(z/Lp) 
Li 


and  the  total  area  can  be  given  by 

A(z)  =  Ax(z)  +  Av(z) 


(36) 


(37) 


At  this  point,  it  is  convenient  to  relate  A^  and  R  with  a  geometrically 
determined  dimensionless  variable,  (3,  by  the  equation 

Ax  =  R2R2  (38) 

Now  by  substituting  equations  (36)  through  (38)  into  equation  (35)  and 
solving,  a  first  order  ordinary  differential  equation  is  obtained  which  relates 
the  radius  of  curvature  to  the  axial  position  along  the  pipe,  or 


dR 

dz 


8nQr  h(z/L  ) 
_ m  p/ 

aL 


v 


l  K'  (A  -  f$2r  )2  K' IS4r  4 
\r  nr  1  m 


(39) 


Since  the  cross-sectional  area,  A,  is  assumed  a  constant,  this  equation  is 
separable  and  therefore  solvable.  The  solution  is 


Q 


max 


0. 16R(K' v)1/2 
8nH(Lp) 


(40) 
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where  H  is  the  integral  of  h  over  the  length  of  the  pipe,  as  presented  by  Cotter 
(1984). 

2.2  Viscous  Limitation 

Vapor  flow  to  the  condenser  relies  on  the  pressure  difference  between 
the  condenser  and  evaporator.  At  low  temperatures,  the  pressure  difference 
across  the  heat  pipe  is  not  great  enough  to  overcome  the  viscous  forces  that 
exist  at  these  conditions.  Dunn  and  Reay  (1983)  found  this  limitation,  called 
the  viscous  limit,  to  be  evaluated  by  the  equation 


%  = 


r2  Lp  P 

V  V  V 


16"vLe» 


(41) 


where  is  the  vapor  pressure  within  the  heat  pipe.  For  most  cases  this 
limit  may  be  neglected. 

2.3  The  Sonic  Limitation 

Sonic  limitations  in  heat  pipes  are  analogous  to  the  sonic  limitations  in 
converging-diverging  nozzles  (Chi,  1976).  In  a  converging-diverging 
nozzle,  the  mass  flow  rate  is  constant  and  the  vapor  velocity  varies 
because  of  the  changing  cross-sectional  area.  In  heat  pipes,  the  reverse 
occurs;  the  area  is  constant  and  the  vapor  velocity  varies  because  of  the 
evaporation  and  condensation  along  the  heat  pipe.  As  in  nozzle  flow, 
decreased  outlet  pressure,  or  in  this  case  condenser  temperature,  results 
in  a  decrease  in  the  evaporator  temperature  until  the  sonic  limitation  is 
reached.  Any  further  increase  in  the  heat  rejection  rate  does  not  reduce 
the  evaporator  temperature  or  the  maximum  heat  transfer  capability,  but 
only  reduces  the  condenser  temperature,  because  of  the  existence  of  choked 
flow. 

The  sonic  limitation  in  heat  pipes  can  be  determined  as 
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q  =  A  P  L 
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^  rR  T  ' 

V  V 


12(y  +  1)J 


(42) 


where  Ty  is  the  mean  vapor  temperature  within  the  heat  pipe. 

2.4  The  Entrainment  Limitation 

Since  the  liquid  and  vapor  flow  in  opposite  directions  in  a  heat  pipe,  at 
high  enough  vapor  velocities,  liquid  droplets  may  be  picked  up  or 
entrained  in  the  vapor  flow.  This  entrainment  results  in  excess  liquid 
accumulation  in  the  condenser  and  hence,  dryout  of  the  wicking  structure. 
Therefore,  it  is  necessary  to  evaluate  when  the  onset  of  entrainment  begins  in  a 
counter-current,  two-phase  flow.  The  most  commonly  quoted  criterion  to 
determine  this  onset  is  when  the  Veber  number,  Ve,  defined  as  the  ratio  of 
the  viscous  shear  force  to  the  force  resulting  from  the  liquid  surface 
tension  or 


We 


2<ch.g>°/ 


v 


a 

is  equal  to  unity.  To  prevent  the  entrainment  of  liquid  droplets  in  the 
vapor  flow,  the  Weber  number  must  therefore  be  less  than  one. 

By  relating  the  vapor  velocity  to  the  heat  transport  capacity 


(43) 


V  = 
v 


(44) 


A  p  L 
v  v 


a  value  for  the  maximum  transport  capacity  based  on  the  entrainment 
limitations  may  be  determined  as 


q  =  A  L 


,1/2 


l  2r 


(45) 


where  r^  is  the  hydraulic  radius  of  the  groove.  A  somewhat  different 
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approach  has  been  proposed  by  Rice  and  Fulford  (1987)  who  developed  an 
expression  that  defines  the  critical  dimensions  for  wicking  structures 
in  order  to  prevent  entrainment. 

2.5  The  Boiling  Limitation 

When  sufficient  heat  is  applied  to  the  evaporator  of  a  heat  pipe, 
boiling  may  occur  and  hence  trap  bubbles  in  the  wick.  This  trapping  of 
bubbles  blocks  the  return  of  fluid  to  the  evaporator  and  causes  dryout  of 
the  wick,  a  phenomenon  referred  to  as  the  boiling  limit.  The  actual 
formulation  of  this  limit  is  based  on  nucleate  boiling  theory.  For  a  heat 
pipe  with  an  axial  evaporator,  this  limit  is  described  by  the  equation 


9k  = 


2rtL  k  T 
e  e  v 


(  2  a 


Lp  ln(r./r  ) 
v  '  1  v 


-  A P 


V  r 


c,m 


(46) 


where  k  is  the  effective  conductivity,  r.  is  the  inner  radius  of  the  heat  pipe, 
rv  is  the  vapor  core  radius  and  rn  is  the  nucleation  radius.  Investigations  by 
Chi  (1976)  have  found  the  nucleation  radius  to  be  within  a  range  of  2.54x10 to 
2.54x10  ^  meters  for  conventional  heat  pipes. 

2.6  Discussion  of  the  Steady-State  Modeling  Results 

To  compare  the  previously  described  limitations,  a  FORTRAN 
computer  code  was  developed.  This  computer  file  calculates  all  the  limits 
previously  mentioned  except  the  boiling  limit.  The  limits  can  be  calculated  as 
a  function  of  the  operating  temperature  or  tilt  angle.  The  working  fluid  can  be 
selected  from  four  working  fluids,  water,  methanol,  ammonia,  or  acetone. 

The  input  dimensions  used  in  the  calculation  of  the  limits  presented  in  the 
following  discussion  are  listed  in  Table  2.  The  calculation  of  Cotter's 
capillary  limit  required  the  several  dimensionless  constants  obtained  from  the 
original  paper.  These  values  are  presented  in  Table  3  for  clarity.  It  is 
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significant  to  note  that  since  one  of  the  assumptions  was  that  the  heat  pipes 
had  been  properly  charged,  the  performance  is  not  affected  by  the  amount  of 
working  fluid  present. 


Table  2  Micro  Heat  Pipe  Dimensions 


length 

condensor  length 
evaporator  length 
w 
a 

vapor  channel  diameter 


0.05  m 
0.0127  m 
0.0127  m 
0.000133  m 
0.5854  rad 
0.0006  m 


Table  3  Dimensionless  Constants 


H(Lp) 

,1 


K' 

K' 


0.5 

0.600 

0.975 

2.044 


As  illustrated  in  Fig.  3,  the  capillary  limit  is  clearly  the  controlling 
limit  for  micro  heat  pipe  operation,  since  both  the  conventional  method  and 
Cotter's  method  of  calculation  of  the  capillary  limit  are  well  below  the  other 
limits.  For  this  reason,  further  discussion  will  focus  on  the  capillary  limit 
only. 

The  development  of  the  two  methods  used  to  calculate  the  capillary  limit 
have  some  significant  differences.  A  few  basic  assumptions  are  shared  between 
the  two  methods.  These  are: 

•  the  heat  pipe  is  properly  charged 

•  the  length  to  diameter  ratio  is  much  greater  than  one 

•  the  cross-sectional  area  is  relatively  constant 

•  the  vapor  and  liquid  are  at  the  same  constant  temperature 
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-  Steady  —  Siate  Capillary  Limit 
•  Capillary  Limit  (Cotter.  1984) 
Sonic  Limit 


OPERATING  TEMPERATURE  (C) 


The  steady-state  model  described  earlier,  makes  no  further  assumptions  while  the 
method  proposed  by  Cotter  (1984),  assumes: 

•  gravity  is  negligible 

•  the  local  axial  heat  transfer  is  specified 

•  the  liquid-vapor  interface  velocity  is  zero 

•  the  vapor  flow  regime  is  both  laminar  and  incompressible 

To  determine  the  significance  of  the  gravitational  forces,  Figs.  4,  5,  and 
6  are  provided.  For  the  case  of  a  horizontal  heat  pipe  with  water  as  the 
working  fluid,  Fig.  4,  the  gravitational  effects  which  appear  in  the  hydrostatic 
pressure  term  as  the  horizontal  axis  are  clearly  negligible.  However,  when  the 
temperature  is  held  constant  and  the  tilt  angle  is  varied,  Fig.  5,  the 
gravitational  effects  become  more  apparent.  Figure  6  illustrates  how  the 
performance  of  a  micro  heat  pipe  as  predicted  by  steady-state  model  diminishes 
when  the  wicking  structure  has  to  work  against  gravity.  Together,  Figs.  5  and  6 
illustrate  that  as  is  the  case  for  conventional  heat  pipes,  for  increasing  tilt, 
the  heat  transport  of  the  heat  pipe  decreases.  This  would  indicate  that  at  a 
tilt  angle  other  than  zero  (horizontal),  gravity  has  a  significant  effect  on  the 
capillary  limit  and  cannot  be  assumed  negligible. 

The  assumption  that  the  local  axial  heat  transport  be  specified,  results  in 
a  parameter,  H(L),  which  must  be  known  a  priori  to  predict  the  maximum  transport 
capacity.  In  this  sense,  it  behaves  as  a  correction  factor  which  can  be  varied 
to  better  simulate  the  experimental  results,  but  is  of  little  value  in  actually 
predicting  the  performance  limitation  prior  to  testing.  Figure  7  illustrates 
the  effect  of  this  parameter  on  the  capillary  limit  as  predicted  by  Cotter 
(1984).  As  shown,  variation  of  this  parameter  results  in  significant  variations 
in  the  predicted  transport  capacity  with  a  value  of  0.25  resulting  in  the 
closest  correlation  between  the  two  techniques. 

Two  additional  assumptions  were  made  in  the  development  of  the  Cotter 


20 


1000 


OPERATING  TEMPERATURE  (C) 
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Comparison  of  the  capillary  limits  as  a  function  of  the  tilt  angle  for 
the  trapezoidal  heat  pipe  with  water  at  25°C. 


model.  First,  is  that  the  liquid-vapor  interface  velocity  is  zero.  Although 
this  introduces  some  error,  it  is  small  and  is  necessary  in  order  to  proceed. 
Second,  is  the  assumption  that  the  vapor  flow  regime  is  both  laminar  and 
incompressible.  In  all  of  the  cases  evaluated,  the  flow  regime  of  the  vapor  was 
found  to  be  in  a  region  that  would  make  this  assumption  both  reasonable  and 
appropriate. 

III.  TRANSIENT  MODEL 

To  understand  the  effects  of  parameters  such  as  the  geometry,  the 
liquid-solid  contact  angle,  and  the  amount  of  working  fluid,  knowledge  of  the 
transient  behavior  is  required.  This  coupled  with  the  small  size,  it  is  very 
difficult  to  measure  the  thermal-hydraulic  parameters  within  a  given  pipe  during 
operation  resulted  in  the  development  of  a  transient  model  capable  of  predicting 
the  behavior  during  start-up  or  during  variations  in  the  thermal  load  in  the 
evaporator  portion  of  the  heat  pipe.  The  steady-state  modeling  predictions  made 
previously  indicate  that  the  boiling,  sonic,  and  entrainment  limits  for  heat 
pipes  of  this  geometry  are  well  above  the  capillary  limits  and  as  a  result,  the 
model  considers  only  the  capillary  limit.  Also,  because  of  the  symmetry  of  the 
heat  pipe,  a  one-dimensional  model  was  used  with  all  of  the  variables  presented 
as  functions  of  x,  the  position  along  the  heat  pipe,  and/or  t,  the  amount  of 
time  since  the  last  change  in  the  evaporator  power  level. 

3.1  Modeling  Assumptions 

To  better  understand  the  capability  and  limitations  of  the  transient  model, 
the  following  is  a  list  of  the  assumptions  made.  We  attempted  to  minimize  these 
in  order  to  improve  the  applicability  of  the  model. 

(1)  Both  the  liquid  and  vapor  flows  are  assumed  to  be  laminar.  Although  the 
model  contains  a  check  on  this,  according  to  computational  results  for  0.2 
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watts  of  power  input  to  the  evaporator,  the  maximum  Reynolds  number 
obtained  for  the  liquid  and  vapor  flow  were  0.22  and  15.4  respectively. 
Both  of  these  values  are  well  below  the  range  of  Re  <  2000. 

(2)  The  evaporation-condensation  rate  is  assumed  to  be  proportional  to  the 
difference  between  the  vapor  pressure  and  the  saturation  pressure  and 
independent  of  the  curvature  of  the  liquid-vapor  interface. 

(3)  The  radius  of  curvature  for  the  liquid  vapor  interface  is  assumed  to  be  a 
function  only  of  the  geometry  of  the  heat  pipe,  the  wetting  angle  and  the 
liquid  cross-sectional  area. 

(4)  The  vapor  within  the  vapor  portion  of  the  heat  pipe  is  assumed  to  behave 
as  an  ideal  gas.  This  assumption  is  reasonable  because  of  the  very  low 
vapor  pressures  involved. 

(5)  The  liquid  distribution  within  each  of  the  four  corners  of  the  heat  pige 
is  assumed  to  be  uniform.  Evaluation  of  the  bond  number  (Bq  =  c^/p^gl 
which  is  the  ratio  of  the  viscous  forces  to  the  gravitational  forces  is 
much  larger  than  unity  for  the  geometry  considered  here.  This  implies 
that  the  viscous  forces  dominate. 

(6)  Computation  of  the  liquid-vapor  interface  is  accomplished  by  treating  this 
interface  as  a  smooth  solid  boundary. 

(7)  The  transfer  of  heat  at  the  liquid  vapor  interface  due  to  conduction  is 
assumed  to  be  negligible. 

3.2  Governing  Equations 

The  high  thermal  conductivity  of  heat  pipes  is  the  result  of  the  evaporation 
and  condensation  process  occurring  within  the  heat  pipe  and  the  high  latent  heat 
associated  with  these  processes.  For  this  reason,  determination  of  the 
evaporation  and  condensation  rate  of  the  liquid  plays  a  key  role  in 
determination  of  the  heat  pipe  behavior.  The  equation  for  the  molecular  flux  of 
evaporation,  j,  given  by  Collier  (1972)  was  employed  in  this  model.  This 
relationship  can  be  expressed  as 
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where,  M  is  the  molecular  weight  of  the  fluid,  R  is  universal  gas  constant,  T  is 
the  absolute  temperature,  Pg  is  the  saturation  pressure,  and  is  vapor 
pressure.  When  Pg  is  greater  than  P^,  j  is  positive  and  the  liquid  evaporates. 
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When  P  is  less  than  P  ,  j  is  negative  and  the  vapor  condenses.  In  the 

S  v 

numerical  model,  the  evaporation-condensation  rate  is  assumed  to  be  proportional 
to  the  liquid-vapor  interface  area  in  each  section  of  the  heat  pipe,  i.e. 


AM  =  j  M  =jv  Ax 
ec  e  J  ec 


(48) 


where  AA  is  the  evaporation-condensation  surface  area  in  one  section,  w  is 
ec  ec 

the  perimeter  of  the  liquid-vapor  interface,  and  Ax  is  the  section  length. 

With  a  given  latent  heat,  L,  the  rate  of  the  heat  released  or  absorbed  in 
any  section  can  be  calculated  by  the  equation 


AQ  =  LAM 
ec  ec 


(49) 


For  the  heat  pipe  to  operate  properly,  the  conditions  existing  between  the 
capillary  pressure,  liquid  pressure  and  the  vapor  pressure  must  be  such  that 
flow  of  the  liquid  back  to  the  evaporator  is  promoted.  This  equation  is  the 
same  for  the  micro  heat  pipe  as  for  larger  heat  pipes  and  can  be  expressed  as 


%1  =  Pm  -  P!  = 


(50) 


where,  is  the  vapor  pressure,  a  is  the  liquid  surface  tension,  and  r  is  the 
radius  of  curvature  of  the  liquid-vapor  interface.  Since  both  and  r  are 
functions  of  x,  the  derivative  of  the  above  equation,  yields 
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dx  r  dx 

For  a  finite  interval,  dx  — >  Ax,  hence,  the  finite  form  of  the  equation  can 
be  expressed  as 
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where 
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(53) 


AP  ,  =  AP, 
vl  fv 


AP  + 
mv 


APfl  +  AP 


ml 


and  the  pressure  changes  are  those  due  to  the  liquid  and  vapor  friction  and  the 
momentum  changes  in  one  section,  Ax,  along  the  heat  pipe. 

Based  on  the  momentum  equation,  the  pressure  drops  due  to  friction  for 
laminar  flow  can  be  expressed  as 

2uAxV 

APf  = - —  (54) 

rh 

where  u  is  the  viscosity  of  the  fluid,  r^  is  the  hydraulic  radius,  and  V  is  the 
velocity  of  the  flow  and  is  equal  to  M/pA. 

The  pressure  drop  due  to  momentum  change  is  given  by 


VAM  +  MAV  +  AM  AV 


AP  =- 
m 


(55) 


where  M  is  the  mass  flowrate  and  A  is  the  flow  cross-sectional  area.  In  this 

model,  AM  equals  AM  for  the  vapor  and  -AM  for  the  liquid. 

ec  ec 

The  energy  equation  used  in  this  model  takes  the  form 


where 


Qin 


out 


+  Q_ 


(56) 


P1^1^P1  ^ 


p  A  C„  AT  dx  + 
\  v  Pw  v 


p  A  C„  AT  dx 
s  s  Ps  s 


(57) 


Cpi’  cpv  anc*  cps  are  t*ie  ^eat  capacities  of  the  liquid,  vapor  and  solid 
boundaries  respectively  and  the  AT's  are  the  corresponding  temperature  changes. 
The  continuity  equation  takes  the  form 


p  A  dx  + 
v  v 


p-^  A^  dx  =  constant  =  fluid  charge 


(58) 


and  total  volume  inside  the  heat  pipe  can  be  defined  and  expressed  as: 
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(59) 


Vol  =  Vol.  +  Vol  and  A  =  A,  +  A 
1  v  1  \> 

3.3  Modeling  Approach 

As  is  the  case  for  larger  heat  pipes,  the  micro  heat  pipe  also  consists  of 
three  basic  sections,  an  evaporator  and  a  condenser,  separated  by  an  adiabatic 
section.  Each  of  these  sections  has  a  different  set  of  boundary  conditions  and 
needs  to  be  treated  accordingly.  For  this  reason,  they  will  be  discussed 
separately  in  the  following  sections. 

3.3.1  Evaporator  Region  The  single  boundary  condition  utilized  in  the 
evaporator  section  is  the  input  heat  density  as  a  function  of  time,  t,  and 
position,  x.  For  a  given  input  heat  density  and  vapor  pressure  at  any  location, 
x,  the  saturation  pressure  can  be  obtained  by  using  the  combination  of  molecular 
flux  equation  (47)  and  the  energy  conservation  equation  (57)  at  a  given 
position.  These  equations  assume  the  following  form 


AQ 

- ££ -  +  p 

(M/2nRT1)1/^  Axu>ec  V 

Op.  P  A.  Ax  AT. 

AQ.  =  — — - - - - - 

At 


(60) 


(61) 


AQ.  =  AQ  +  AQ.  (62) 

in  ec  1 

Since  Tj  is  a  function  of  Pg,  equations  (60  through  62)  can  be  correlated  to 

each  other.  An  iteration  method  with  relaxation  was  used  to  solve  the  system  of 

equations  to  find  the  value  of  T,,  Q  and  Q.  .  The  difference  between  the 

l  ec  i 

boundary  and  liquid  temperatures,  AT^,  is  proportional  to  the  input  heat 
density.  The  boundary  temperature  can  therefore  be  obtained  by  adding  AT^  to 
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3.3.2  Adiabatic  Region.  In  the  adiabatic  section,  there  is  by  definition  no 
heat  into  or  out  of  the  heat  pipe.  The  heat  absorbed  by  evaporation  or  released 
by  condensation  is  equal  to  the  heat  consumed  in  heating  up  or  cooling  down  the 
liquid  and  the  heat  pipe  shell.  The  governing  equations  for  this  section 
therefore  become 


Q  =  Leo  Ax 
ec  ec 


M  U/2 


l  2nRT1  ) 


(Ps  -  Pv) 


Q1  =  CpipiAiATiAx 


Q  =  Cn  p  A  AT,  Ax 
s  Ps  s  s  b 


<2  +  Qi  +  Cl  *  0 

ec  i  s 


(63) 

(64) 

(65) 

(66) 


Because  the  magnitude  of  Qgc  in  the  adiabatic  section  is  much  less  than  that 
in  the  other  two  sections,  the  heat  transferred  between  liquid  and  the  heat  pipe 
shell  is  very  small,  and  hence  the  temperature  drop  can  be  considered 
negligible.  By  assuming  T^  =  T^,  the  equations  given  above  can  be  solved  by 
iteration,  to  obtain  Q  ,  Q^j  Qg  and  T^. 

3.3.3  Condenser  Region.  In  the  condenser  region,  the  boundary  temperature 
of  the  heat  pipe  is  assumed  to  be  fixed.  Therefore,  the  governing  equations  are 
similar  to  those  used  in  the  adiabatic  section 


Q  =  Lw  Ax 
ec  ec 


l  2  JlRT1  J 


1/2 


(Ps  -  Pv) 


Qj  —  Cp^p^A^AT^Ax 


(67) 

(68) 


out 


=  Q  + 
ec 


(69) 


30 


where  AT^  is  proportional  to  the  heat  output.  For  a  given  T^,  these  equations 
are  solved  for  Qec»  and  using  an  iterative  technique. 

3.4  General  Assembly 

Although  the  heat  pipe  is  small,  many  variables  are  involved  in  the 
determination  of  the  transient  behavior.  Some  of  the  most  important  ones  are 
listed  below: 


—  Liquid  cross-section  area 
—  Vapor  cross-section  area 
M  —  Evaporation-Condensation  Rate 
M-^  —  Liquid  Hass  Flowrate 

—  Vapor  Mass  Flowrate 
P  —  Vapor  Pressure 

P-^  —  Liquid  Pressure 

Q  —  Heat  Input  and  Output 

r  —  Radius  of  Curvature  on  Liquid-Vapor  Interface 

—  Boundary  Temperature 
T^  —  Liquid  Temperature 
T^  —  Vapor  Temperature 

All  these  parameters  are  functions  of  both  x  and  t.  In  addition,  most  of 
them  are  interrelated  and  therefore  solutions  cannot  be  obtained  explicitly. 

This  serves  to  significantly  increase  the  complexity  of  the  model.  A  double 
iteration  technique  is  employed  to  minimize  the  error.  One  is  false  position, 
which  is  used  to  find  the  root  of  the  implicit  function.  The  other  is  a 
relaxation  technique  used  to  find  solutions  for  the  implicit  array. 

Since  the  evaporation-condensation  rate  is  very  sensitive  to  changes  in  the 
vapor  pressure,  a  special  treatment  is  applied  to  obtain  a  stable  solution. 
First,  the  evaporation  rate  is  assumed  to  be  equal  to  the  condensation  rate  and 
the  vapor  pressure  is  calculated.  Second,  the  vapor  pressure  is  used  to  compute 
the  change  in  the  vapor  mass  and  hence,  the  difference  between  evaporation  and 
condensation.  The  original  assumption  of  equal  evaporation  and  condensation 
rates  are  adjusted  and  updated.  Preliminary  tests  indicate  that  the  solution  is 
stable  regardless  of  the  time  step  size  used.  In  comparison,  prior  to  using  the 
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predictor-corrector  treatment,  the  solution  vibrated  with  increasing  magnitude 
for  time  step  sizes  below  0.0001  second. 


The  general  computing  procedures  are  listed  below: 


1.  Initialize  the  variables  and  arrays. 

2.  Set  the  time  equal  to  the  first  time  step. 

3.  For  a  given  boundary  condition,  compute  the  boundary,  liquid  and  vapor 
temperatures  and  the  evaporation-condensation  rate  for  each  section  by 
assuming  that  the  total  evaporation  equals  the  total  condensation. 

4.  Compute  the  pressure  drops  and  flowrates  for  each  section  by  iteration. 

5.  Compute  the  radii  of  curvature  for  the  liquid-vapor  interface. 

6.  Compute  the  vapor  pressure  along  with  the  boundary,  liquid  and  vapor 
temperatures  and  the  evaporation-condensation  rate  for  each  section  by 
i teration. 

7.  Make  a  correction  on  the  liquid-vapor  mass  change  using  the  predicted 
vapor  pressure. 

8.  Repeat  steps  4  through  6  with  the  vapor  pressure  correction. 

9.  Calculate  the  liquid  and  vapor  flowrate  once  again  with  the  corrected 
vapor  pressure.  Calculate  the  liquid  pressure. 

10.  Renew  all  of  the  variables  and  arrays  and  print  out  the  results. 

11.  If  the  time  is  less  than  the  given  limit,  increase  on  time  step  and  go 
back  to  step  2. 

3.5  Discussion  of  the  Transient  Modeling  Results 

3.5.1  Geometry.  Because  both  the  steady-state  and  transient  performance  of 
the  heat  pipe  are  strongly  dependent  upon  the  physical  geometry,  it  is  necessary 
to  accurately  describe  this  geometry.  The  actual  shape,  shown  in  Fig.  8,  the 
shape  can  be  described  as  follows 


d  ,  n 

a  =  -  cos  -  (70) 

2  4 
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d. 

h  =  a  -  —  (71) 

2 


r  ■=  - 

S  2h 

(72) 

0  =  sin-*  (a/rs) 

(73) 

Xg  =  -rg  sin  (n/4  -  n) 

(74) 

yQ  =  rg  cos  (n/4  -  n) 

(75) 

=  ^0  -  ^  rs2  -  (x  -  x0)2 

(76) 

2 

Aarc  =  rs  ^  “  sin0  cos0) 

(77) 

A  =  (d,  sin  n/4)^  -  4A 

d  '  arc 

(78) 

A  =  d  2  -  A 
s  o 

(79) 

dy  x  -  x0 

dX  J  2  ■  ,2 
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a  =  tan  *  (dy/dx) 
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The  dimensions  obtained  for  the  pipe  tested  were  dQ=  1.0  mm,  d^=  0.6  mm,  and 

d,=  1.0  mm. 
d 

Once  the  geometry  was  known,  it  was  necessary  to  determine  the  contact 
angle,  9,  which  was  computed  as: 


0  =  cos 


a  -  a  , 
sv  si 


3lg 


(94) 


Although  some  information  is  available  on  the  surface  tension  between  the 
solid  and  the  liquid,  asp  only  limited  information  on  the  other  two  values,  the 
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surface  tension  between  the  solid  and  the  vapor,  ctsv,  and  the  liquid  and  the 
vapor,  c-^v,  was  available.  For  this  reason,  an  experimental  investigation  was 
conducted  to  determine  the  contact  angle  between  copper  and  water,  and  also 
between  silver  and  water. 

An  acrylic  chamber  with  internal  dimensions  of  approximately  2.54  cm  by  2.54 

cm  by  8.0  cm  was  constructed.  The  front  and  back  walls  of  the  chamber  were 

transparent  while  the  left  and  right  sides  were  made  of  copper  and  silver 

respectively.  The  chamber  was  filled  through  a  fill  port  located  in  the  top, 

3 

with  approximately  10  cm  of  distilled  water.  The  chamber  was  evacuated  for 
approximately  5  minutes  to  eliminate  noncondensible  gases.  The  fill  port  was 
then  closed,  sealing  the  chamber.  By  tilting  the  chamber  until  the  liquid  at 
the  interface  formed  a  flat  surface  as  shown  in  Fig.  9,  the  contact  angle  could 
be  measured.  The  results  of  these  experimental  measurements  indicated  contact 
angles  of  55°  and  45°  for  copper  and  silver  respectively. 


Fig.  9  Photograph  of  the  contact  angle  experimental  device. 
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The  preliminary  results  from  the  numerical  model  indicated  that  this 
parameter,  the  contact  angle,  was  of  significant  importance  in  determining  the 
liquid  pressure  drop,  the  maximum  capillary  pumping  pressure,  the  maximum  heat 
transport  capacity,  and  the  time  required  to  reach  steady-state. 

Figure  10  compares  the  liquid  cross-sectional  area  with  the  liquid-vapor 

interfacial  pressure  drop  for  wetting  angles  of  0°  and  45°.  As  illustrated,  the 

2 

maximum  liquid  cross-sectional  area  for  a  contact  angle  of  45°  is  0.218  mm  , 

2 

much  larger  than  the  0.116  mm  for  a  contact  angle  of  0°.  As  a  result,  the 
difference  between  the  minimum  interface  curvature  in  the  evaporator  and  the 
maximum  interface  curvature  in  the  condenser  is  much  larger  for  the  45°  contact 
angle  than  for  the  case  of  the  0°  contact  angle  for  the  same  heat  input,  as 
illustrated  in  Fig.  11.  Near  dryout,  the  frictional  pressure  drop  in  the 
evaporator  is  much  larger  for  the  45°  contact  angle  than  for  the  0°  case  as 
illustrated  in  Fig.  12.  The  resulting  effects  of  variations  in  the  contact 
angle  are  large  differences  in  the  liquid  pressure  drop  and  the  capillary 
pumping  pressure. 

The  relationships  described  above  do  not  necessarily  mean  that  large  contact 
angles  will  result  in  a  larger  maximum  heat  transport  capacity,  since  the  heat 
transport  capacity  is  also  dependent  upon  the  liquid  charge  within  the  heat 
pipe.  Using  the  fill  charge  of  the  actual  test  pipes,  and  a  31°C  condenser  bath 
temperature,  a  maximum  heat  transport  capacity  of  0.11  watt  for  a  45°  contact 
angle  resulting  from  evaporator  dryout  was  predicted  as  compared  to  a  maximum 
value  of  0.12  for  a  0°  contact  angle  resulting  from  condenser  flooding.  In 
addition  to  affecting  the  pressure  drops,  the  transient  analysis  indicates  that 
the  time  required  for  the  heat  pipe  to  reach  steady-state  is  longer  for  larger 
contact  angles.  Once  a  heat  input  of  0.1  watt  has  been  reached,  a  total  of  80 
seconds  is  required  for  the  heat  pipe  to  reach  steady-state  for  contact  angles 
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terfacial  pressure  drop  as  a  function  of  the  liquid  cross-sectional  area. 


45  Degree  Contact  Angle 
0  Degree  Contact  Angle 


Fig.  11  Steady-state  liquid  cross-sectional  area  as  predicted  by  the  transient  model. 


of  45°,  while  only  30  seconds  are  required  if  the  contact  angle  is  0°. 

Because  of  the  geometrical  shape  of  the  heat  pipe,  a  phenomena  referred  to 
as  an  inverted  miniscus  occurs  at  contact  angles  greater  than  45°.  This 
phenomena  can  best  be  described  as  a  change  in  the  shape  of  the  liquid-vapor 
interface  from  one  that  is  convex  to  one  which  is  concave.  From  the  above 
analysis,  it  is  apparent  that  the  contact  angle  is  one  of  the  key  parameters  in 
the  determination  at  the  optimum  charge  and  the  heat  transport  capacity  of  the 
heat  pipe. 

The  steady-state  temperature  profile  as  predicted  by  the  transient  model  is 
illustrated  in  Fig.  13  for  an  input  power  of  0.1  watt  and  a  contact  angle  of 
45°.  As  shown,  the  liquid,  vapor,  and  boundary  temperatures  are  all  the  same  at 
the  left  end  of  the  pipe.  This  point  represents  the  highest  temperature  within 
the  pipe.  The  vapor  temperature  decreases  slightly  in  the  evaporator  region  and 
then  remains  essentially  constant  throughout  the  adiabatic  and  condenser 
regions.  The  boundary  and  liquid  temperatures  decrease  somewhat  more  rapidly  in 
the  evaporator  region  and  continue  to  decrease  in  both  the  adiabatic  and 
condenser  regions.  The  two  large  gradients  are  apparent  at  the  exit  to  the 
evaporator  and  the  entrance  to  the  condenser  and  are  the  result  of  longitudinal 
step  size  used  in  the  model.  Because  the  wall  of  the  heat  pipe  is  very  thin, 
the  temperature  difference  between  the  wall  and  the  liquid  temperatures  are  very 
small. 

The  transient  model  indicates  that  the  vapor  temperature  variations  are 
proportional  to  the  heat  input,  and  that  without  dryout  or  flooding,  the 
conductivity  of  the  heat  pipe  is  close  to  a  constant  value  with  respect  to  time. 
In  contrast  to  the  vapor  temperature,  a  fairly  long  period  of  time  is  required 
for  the  other  thermal  hydraulic  parameters  to  reach  their  steady-state  values. 
The  time  period  for  the  case  mentioned  previously  ranges  from  30  seconds  to  80 
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Fig.  12  Steady-state  liquid  and  vapor  pressure  as  a  function  of  position  as  predicted  by  the  transient  model 
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y-state  temperature  distribution  as  predicted  by  the  transient  model. 


seconds  after  full  heat  input,  depending  on  the  contact  angle.  The  smaller  the 
contact  angle  and  the  higher  the  temperature,  the  shorter  the  time  required. 

Figures  14  and  15  illustrate  the  changes  in  the  liquid  mass  flow  rate  and 
liquid  cross-sectional  area  as  a  function  of  time.  It  is  significant  to  note 
that  during  the  time  immediately  after  the  application  of  the  heat  to  the 
evaporator,  the  liquid  and  vapor  flow  in  the  evaporator  are  both  in  the  same 
direction.  This  reverse  liquid  flow  is  the  result  of  an  imbalance  in  the  total 
pressure  drop.  The  evaporation  rate  does  not  provide  a  change  in  the 
liquid-vapor  interfacial  curvature,  adequate  to  compensate  for  the  total 
pressure  drop  resulting  from  the  temperature  change.  Once  the  heat  input 
reaches  full  load  the  reverse  liquid  flow  disappears  and  the  liquid  mass  flow 
rate  gradually  increases  until  it  reaches  a  steady-state  condition  in  which  the 
liquid  mass  flow  rate  is  equal  to  the  vapor  mass  flow  rate  in  any  given  section. 
This  is  more  apparent  in  Fig.  16  which  illustrates  the  pressure  as  a  function  of 
axial  position  for  both  10  seconds  after  start-up  and  for  steady-state.  As 
illustrated,  at  10  seconds  after  start-up,  the  pressure  of  both  the  liquid  and 
vapor  are  higher  in  the  evaporator  and  gradually  decrease  with  position, 
promoting  flow  away  from  the  evaporator,  while  at  steady-state,  the  liquid 
pressure  is  much  lower  in  the  evaporator  and  increases  dramatically  with 
position,  promoting  flow  back,  to  the  evaporator. 

Finally,  it  is  important  to  note  that  the  pressure  drops  due  to  friction  are 
much  larger  than  for  either  the  pressure  drops  due  to  the  vapor  momentum  change 
(two  orders  of  magnitude)  or  the  pressure  drop  due  to  the  liquid  momentum  change 
(four  orders  of  magnitude). 
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Liquid  mass  flow  rate  as  a  function  of  time  and  the  axial  position 
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Liquid  cross-sectional  area  as  a  function  of  time  and  the  axial  position 
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Liquid  and  vapor  pressure  as  a  function  of  time  and  the  axial  position 


IV.  EXPERIMENTAL  PROGRAM 

4.1  Description  of  Test  Article 

The  test  articles  used  in  this  investigation  were  obtained  from  ITOH 
Research  and  Development  Corporation  of  Japan.  The  physical  size  and  internal 
dimensions  are  illustrated  in  Figs.  1  and  2  respectively.  As  shown,  the  outer 
dimensions  are  1  mm  by  1  mm  and  the  length  is  57  mm.  The  case  material  was 
either  copper  or  silver  and  the  working  fluid  used  was  ultrapure  water.  Two  of 
the  test  pipes,  one  silver  and  one  copper,  were  charged  with  0.0032  gram  of 
fluid  prior  to  sealing.  Two  additional  uncharged  test  pipes  were  also  obtained 
to  establish  a  baseline  value  and  determine  the  improvement  in  the  heat 
transport  capacity  due  to  the  vaporization  and  condensation. 

4.2  Description  of  Experimental  Apparatus 

To  conduct  the  experimental  portion  of  the  investigation,  a  test  facility 
was  constructed  in  the  Conduction  Heat  Transfer  Laboratory  at  Texas  A&M 
University.  This  test  facility  illustrated  in  Fig.  17,  consisted  of  five 


Fig.  17  Photograph  of  the  micro  heat  pipe  test  facility  with  the  IR  camera. 


Fig.  18  Mic«.o  heat  pipe  test  stand. 

subsystems:  the  test  stand,  a  variable  heat  source,  a  constant  temperature 

cooling  source,  a  data  acquisition  system,  and  a  vacuum  system.  The  test  stand, 
shown  in  Fig.  18,  was  constructed  from  Plexiglass  and  allowed  the  micro  heat 
pipe  to  rotate  through  an  angle  of  360°.  This  allowed  testing  at  tilt  angles 
which  both  helped  and  hindered  the  return  of  liquid  back  to  the  evaporator.  The 
test  stand  also  supported  an  aluminum  cooling  chamber  which  sunounded  the 
condenser  region  of  the  heat  pipe.  A  Forma  Scientific  circulating  bath  provided 
a  constant  temperature  fluid  to  the  cooling  chamber  where  it  flowed  over  the 
condenser  portion  of  the  heat  pipe.  The  bath  was  controlled  in  such  a  manner 
that  a  coolant,  consisting  of  an  ethyl  glycol  solution  of  constant  temperature 
and  flow  rate  was  continually  provided. 
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The  evaporator  portion  of  the  heat  pipe  was  surrounded  by  a  piece  of  solid 
copper  0.635  cm  in  diameter  and  1.27  cm  in  length.  This  copper  heat  spreader 
was  in  turn  encased  in  an  electric  resistance  heater  manufactured  by  Watlow 
Manufacturing  Company.  The  power  was  controlled  by  a  regulated  power  supply  and 
monitored  by  measuring  the  voltage  and  current  to  the  resistance  heater. 

Two  separate  techniques  were  used  to  monitor  the  temperature  distribution  of 
the  heat  pipe.  The  first  of  these  consisted  of  a  series  of  Chromel-Alumel 
K-type  thermocouples  attached  to  the  outer  wall  of  the  heat  pipe.  Three 
thermocouples  were  located  at  equidistant  positions  on  the  evaporator  portion  of 
the  heat  pipe:  one  at  the  end  of  the  pipe,  one  at  the  end  of  the  evaporator 
portion,  and  one  midway  between.  Two  additional  thermocouples  were  located  in 
the  condenser  region:  one  on  the  outer  surface  of  the  heat  pipe  at  the  entrance 
to  the  condenser,  and  the  other  located  at  the  entrance  to  the  condenser  bath. 

In  addition,  one  thermocouple  was  located  on  the  outer  surface  of  the  heater  for 
use  in  calibration.  The  thermocouples  were  monitored  using  an  HP  3497A  data 
acquisition  system  in  conjunction  with  an  HP  87XM  Personal  Computer  which 
monitored  not  only  the  temperatures  but  the  voltage  and  current  measurements  as 
made  by  two  HP  3478A  multimeters. 

The  second  technique  used  to  monitor  the  temperature  distribution  in  the 
heat  pipe  was  an  infrared  thermal  measurement  system.  A  Hughes  Probe  Eye  TVS 
Model  3000  Infrared  Camera  in  conjunction  with  an  RGB  Monitor  and  a  Real  Time 
Recorder  were  used  to  observe  the  test  pipe  during  both  start-up  and 
steady-state  operation.  Using  this  system  the  temperature  distribution  along 
the  heat  pipe  could  be  measured  with  an  accuracy  of  ±0.05°C. 

To  reduce  any  heat  losses  through  convection,  all  of  the  tests  were 

_2 

conducted  in  a  vacuum  of  >  10  torr.  The  entire  test  facility  was  encased  in  a 
stainless  steel  vacuum  chamber  equipped  with  a  single  crystal  silicon  window 
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which  allowed  observation  with  the  thermal  measurement  system.  A  mechanical 
roughing  pump  in  series  with  an  oil  diffusion  pump  were  used  to  obtain  the 
required  vacuum. 

4.3  Experimental  Procedure 

To  fully  understand  the  operational  characteristics  and  performance 
limitations  of  the  heat  pipe,  we  developed  an  experimental  test  plan  comprised 
of  two  tasks: 

i.  Testing  of  a  charged  copper  pipe  to  determine  the  time  required  for  the 
heat  pipe  to  reach  steady-state. 

ii.  Testing  of  several  charged  pipes,  one  copper  and  one  silver,  to  obtain 
data  with  which  to  determine  the  accuracy  of  the  steady-state  and 
transient  model  (at  steady-state);  to  determine  the  maximum  steady-state 
transport  capacity;  and  to  determine  the  effect  of  variations  in  the  type 
of  case  material  used  and  the  fluid  contact  angle. 

To  accomplish  these  tasks,  an  experimental  procedure  was  developed  and  used 
in  all  of  the  measurements  made.  This  procedure  can  be  described  as  follows: 
once  the  heat  pipe  had  been  instrumented  with  the  thermocouples,  fitted  with  the 
evaporator  heat  sink  and  condenser,  and  mounted  on  the  fixture,  the  heat  pipe, 
evaporator,  and  condenser  were  all  painted  with  a  flat  black  paint  of  known 
emissivity  to  enhance  the  IR  resolution.  The  test  fixture  was  placed  in  the 
vacuum  chamber  and  the  constant  temperature  bath  adjusted  to  the  desired 
temperature.  The  bell  jar  was  lowered  and  the  vacuum  pump  was  turned  on.  The 
heat  pipe  and  the  test  facility  were  allowed  to  outgas  for  a  minimum  of  12 
hours . 

Once  the  desired  vacuum  had  been  obtained,  the  input  heat  flux  was  increased 
in  small  increments  until  the  desired  operating  temperature  had  been  reached. 
Adequate  time  was  allowed  between  each  incremental  increase  to  allow  the  heat 
pipe  to  stabilize.  When  the  desired  operating  temperature  had  been  reached,  the 
heat  addition  continued  until  steady-state  was  reached  (approximately  20  minutes 
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as  illustrated  in  Fig.  19)  at  which  time  the  experimental  data,  including  the 
time,  date,  power  level,  thermocouple  temperatures,  and  coolant  bath  temperature 
were  recorded.  For  the  transient  data  shown  in  Fig.  19  and  for  several 
steady-state  tests,  the  thermal  history  of  the  pipe,  as  measured  by  the  IR 
thermography  system,  was  also  recorded. 

To  obtain  the  data  for  the  next  successive  operating  temperature,  the  heat 
flux  was  incremented  and  the  temperature  of  the  coolant  bath  adjusted  to 
maintain  a  constant  temperature  at  the  entrance  to  the  heat  pipe  condenser.  This 
procedure  was  followed  until  the  conductance  of  the  heat  pipe  approached  a 
constant  value.  At  this  time,  the  temperature  of  the  coolant  bath  was 
incremented,  and  the  process  was  repeated  for  a  new  operating  temperature. 

V.  RESULTS 

Figure  19  compares  the  evaporator  temperature  with  the  elapsed  the  time  for 
two  different  power  increments.  As  shown,  the  heat  pipe  reaches  steady-state, 
for  both  power  levels  quite  rapidly,  with  the  time  required  increasing  with  the 
size  of  the  power  level  increment.  Using  these  data,  a  time  of  1,200  seconds 
(20  minutes)  was  established  as  the  minimum  acceptable  time  between  readings. 
This  ensured  that  the  heat  pipe  would  always  be  operating  at  steady-state  when 
data  were  recorded  and  helped  to  accelerate  the  testing  process  by  eliminating 
the  need  to  continually  monitor  the  heat  pipe  for  temperature  changes. 

Figure  20  illustrates  the  results  obtained  for  a  copper  heat  pipe  of  the 
dimensions  shown  in  Figs.  1  and  2,  operating  in  a  horizontal  orientation  and 
charged  with  0.0032  gm  of  ultrapure  water.  The  heat  pipe  was  evaluated  at  a 
total  of  six  different  operating  temperatures:  30.8°C,  39.5°C,  43.8°C,  51.6°C, 
58.9°C,  and  67.6°C.  The  conductance  value  shown  as  the  vertical  axis  was 
computed  based  upon  the  input  power  to  the  pipe  as  measured  by  the  two 
multimeters.  It  is  important  to  note  the  magnitude  of  the  parameters  compared 
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Experimental  determination  of  the  time  required  to  reach  steady-state  at  different  power  levels 


in  this  figure.  Although  the  variations  in  the  thermal  conductance,  shown  on 
the  ordinate,  appear  very  large,  in  reality  they  are  quite  small  compared  to  the 
changes  in  operating  temperature. 

As  illustrated  in  Fig.  20,  the  thermal  conductance  of  the  pipe  decreased 
rapidly  with  respect  to  the  evaporator  temperature  and  asymptotically  approached 
a  constant  value  somewhat  above  the  theoretical  value  for  an  uncharged  pipe. 

The  conductance  value  then  decreased  only  slightly  with  increases  in  the  input 
power.  With  the  exception  of  the  constant  difference  between  the  ultimate 
conductance  value  of  the  test  pipe  and  the  theoretical  conductance  of  the 
uncharged  pipe,  this  behavior  is  precisely  what  would  be  expected  for  a  heat 
pipe  in  which  the  evaporator  is  slowly  drying  out. 

The  constant  difference  between  the  ultimate  conductance  value  of  the  test 
pipe  and  the  theoretical  conductance  of  the  uncharged  pipe,  can  be  explained  as 
follows:  increases  in  the  input  power  initially  cause  the  liquid  meniscus  to 
recede  into  the  liquid  channels  located  in  the  corners  of  the  heat  pipe.  This 
results  in  a  decrease  in  the  radius  of  curvature  of  the  liquid-vapor  meniscus,  a 
corresponding  decrease  in  the  cross-sectional  area  of  the  liquid  and  hence,  a 
slight  increase  in  the  evaporator  temperature  (all  of  which  are  apparent  from 
the  results  of  the  transient  model,  Figs.  11,  15,  and  16).  The  recession  of  the 
meniscus  continues  until  dryout  of  the  evaporator  occurs.  Once  the  evaporator 
begins  to  dry  out,  the  conductance  of  the  test  pipe  rapidly  approaches  the 
conductance  of  an  uncharged  pipe  as  shown  in  Fig.  20.  Although  the  tendency  is 
the  same  for  all  six  operating  temperatures,  the  approach  toward  the  uncharged 
pipe  value  is  gradual.  This  gradual  approach  is  the  result  of  an  increase  in 
the  effective  size  of  the  evaporator.  As  described  previously,  the  heat  is 
supplied  to  the  heat  pipe  by  an  electrical  resistance  heater  wrapped  around  1.27 
cm  of  the  heat  pipe.  Prior  to  dryout,  all  of  the  evaporation  occurs  over  this 
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Fig.  20  Measured  thermal  conductance  of  a  trapezoidal  micro  heat  pipe  as 
function  of  the  evaporator  temperature  (copper,  0.0032  gm  charge 


region.  Once  dryout  begins,  the  temperature  of  the  heat  pipe  case  in  the 
adiabatic  section  increases,  and  a  small  part  of  what  was  formerly  the  adiabatic 

section  begins  to  act  as  an  evaporator.  As  the  power  continues  to  increase,  the 

dryout  in  the  evaporator  spreads  and  the  portion  of  the  adiabatic  section  which 
behaves  as  an  evaporator  also  expands.  In  the  extreme  case  where  complete 
dryout  of  the  evaporator  section  occurs,  a  large  portion  of  the  adiabatic 
section  may  behave  as  an  evaporator  and  may  even  begin  to  dry  out  at  the  end 
farthest  from  the  condenser. 

Figure  21  illustrates  the  results  obtained  for  a  silver  heat  pipe  of  the 
same  dimensions,  charged  with  0.0032  gm  of  ultrapure  water  and  operating  under 
similar  conditions.  In  this  case,  however,  the  operating  temperatures  at  which 
the  heat  pipe  was  evaluated  were  31.8°C,  32.3°C,  42.0°C,  52.2°C,  and  62.8°C.  As 

shown,  the  trends  and  tendencies  are  all  similar  to  those  of  the  copper  pipe 

tested. 

5.1  Comparison  of  Analytical  and  Experimental  Results 

To  compare  the  results  of  the  experimental  investigation  with  those 
predicted  by  the  analytical  model,  it  was  necessary  to  develop  a  method  for 
determining  the  point  at  which  dryout  begins.  The  experimental  data  presented 
in  Figs.  20  and  21,  clearly  display  a  change  in  the  operation  of  the  heat  pipes 
as  the  operating  temperature  and  power  levels  were  increased.  After  reviewing 
the  experimental  data,  we  attempted  to  identify  and  better  quantify  the  onset  of 
the  dryout  condition.  Because  of  the  expansion  of  the  evaporator  portion  into 
the  adiabatic  section  of  the  heat  pipe  described  previously,  this  proved  to  be 
quite  difficult. 

Several  strategies  were  employed,  including  development  of  a  mathematical 
method  for  determining  the  maximum  curvature  for  each  operating  temperature. 

The  final  approach  used  was  to  select  a  region  on  both  sides  of  the  point  of 
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maximum  curvature,  that  is  a  range  around  the  point  at  which  the  most  rapid 
change  of  the  conductance  occurs  with  respect  to  evaporator  temperature.  This 
range  was  assumed  to  be  the  region  over  which  dryout  occurred. 

The  initial  point,  which  will  be  referred  to  as  the  onset  of  dryout  was 
identified  as  that  point  where  the  slope  of  the  curve  was  60  degrees.  The 
termination  of  dryout  was  assumed  to  occur  at  that  point  where  slope  of  the 
curve  was  30  degrees.  Although  this  would  appear  to  be  a  large  range,  in 
reality  it  is  probably  reasonably  close  to  the  actual  behavior  since  dryout  is  a 
gradual  process  and  unlike  boiling,  develops  quite  slowly. 

Figures  22  and  23  illustrate  the  results  of  this  process  and  compare  the 
experimentally  determined  maximum  heat  transfer  capacity  with  those  predicted  by 
the  steady-state  and  transient  models  for  the  copper  and  silver  test  pipes 
respectively.  When  the  the  condenser  temperature  was  increased  to  temperatures 
above  ambient,  a  small  amount  of  power  was  required  to  maintain  a  uniform 
temperature  along  the  pipe.  To  determine  the  maximum  heat  transfer  capacity, 
shown  on  the  vertical  axis,  this  initial  power  was  subtracted  from  the  actual 
power  input  as  measured  by  the  multimeters. 

As  shown  in  Figs.  22  and  23,  when  the  steady-state  experimental  results  were 
compared  with  the  analytical  model,  they  verified  the  modeling  predictions  for 
both  the  steady-state  model  and  the  transient  model  (at  steady-state).  The 
steady-state  model  overpredicted  the  experimentally  determined  maximum  heat 
transport  capacity  by  approximately  30  percent  for  both  the  copper  and  silver 
test  pipes.  The  transient  model  (at  steady-state),  however,  proved  to  be 
significantly  more  accurate  with  the  difference  between  the  measured  and 
predicted  values  averaging  less  than  10  percent  for  the  copper  test  pipe  and 
5  percent  for  the  silver  test  pipe.  If  a  value  for  H(L)  is  used  in  the  model 
developed  by  Cotter  (1984),  it  predicts  the  onset  of  dryout  reasonably  well, 


57 


o 

o 


o 

CD 


o 

00 


o 

h- 


o 

cO 


O 

if) 


O 


o 

n 


CO  cO 

O  O 


't  K)  CN  <—  O 

o  o  o  6  6 

(m)  U3JSNVdi  1V3H 


(U 

CP 

u 


03 

r 
o  o 
u 

u  e 


•ri  0") 

s 


cn 


05  O 
T3  O 
■H  • 


o  o 

N 

<D  - 

a  m 

03  <1> 


03 

0  (fl 
03 

a 

O  <D 

w  a 

■H  *H 

u  a 

03 

a 

G  03 
O  0) 
U  £ 


CN 

CN 


CT> 

■rl 


58 


OPERATING  TEMPERATURE  (C) 


CD  if)  ^  rO  OJ  <—  O 

6  6  o  o  o  6  o 


(m)  aajSNvyi  ivbh 


59 


OPERATING  TEMPERATURE.  (C) 


although  no  rationale  can  be  made  for  using  0.5  as  opposed  to  1.00  or  0.25, 
which  would  significantly  under  and  overestimate  the  values  respectively. 

In  addition  to  accurately  predicting  the  maximum  heat  transport  capacity, 
the  results  indicated  that  reverse  liquid  flow  occurs  in  the  liquid  channels 
during  start-up  and  that  the  wetting  angle  is  the  single  most  important  factor 
contributing  to  the  transport  capacity  and  behavior  of  these  micro  heat  pipes. 
5.2  Summary 

The  results  of  the  modeling  effort  indicated  that  unlike  traditional  heat 
pipes,  reverse  liquid  flow  occurs  in  the  liquid  channels  of  the  trapezoidal 
micro  heat  pipes  studied  here  during  start-up,  and  that  the  wetting  angle  is  the 
single  most  important  factor  contributing  to  the  transport  capacity  and  behavior 
of  these  types  of  heat  pipes. 

The  steady-state  experimental  results  verified  the  operation  of  the 
trapezoidal  micro  heat  pipe  and  when  compared  with  the  analytical  model  verified 
the  trends  predicted  by  both  the  steady-state  model  and  the  transient  model  (at 
steady-state).  Although  additional  work  is  required  to  refine  the  models, 
either  the  steady-state  or  the  transient  model  (at  steady-state)  could  be  used 
to  accurately  assess  the  effects  of  variations  in  the  physical  geometry,  the 
amount  or  type  of  working  fluid,  and  the  power  levels,  all  at  steady-state.  The 
Cotter  model  (Cotter,  1984)  is  of  limited  value  unless  the  value  for  the 
parameter  H(L)  is  known  in  advance. 

Since  variations  in  the  way  dryout  is  defined  can  significantly  increase  or 
decrease  the  measured  values,  the  discrepancy  between  the  experimentally 
measured  values  and  the  predicted  values  may  in  part  be  due  to  the  technique 
used  to  define  dryout.  As  a  result,  the  percieved  accuracy  of  the  models  is 
strongly  dependent  upon  how  dryout  is  defined.  To  resolve  this  problem,  it  is 
necessary  to  better  understand  the  dryout  phenomena  and  better  define  when  it 
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begins  and  how  it  proceeds. 


Although  the  two  analytical  models  developed  have  been  shown  to  predict  the 
steady-state  performance  limitations  and  operational  characteristics  with  a 
reasonable  degree  of  accuracy,  no  experimental  data  on  the  transient  operational 
characteristics  were  obtained.  In  addition,  several  factors  which  significantly 
affect  both  the  steady-state  and  transient  behavior  require  continued  study. 

First  among  these  is  a  determination  of  the  effect  of  variations  in  the 
cross-sectional  shape  of  the  device;  second  is  the  effect  of  variations  in  the 
amount  of  working  fluid;  and  third,  is  a  determination  of  just  how  small  these 
devices  can  be  and  still  continue  to  function. 

We  recommend  that  both  the  experimental  and  analytical  investigations  be 
continued,  to  address  these  and  other  issues.  The  following  tasks  provide  a 
logical  continuation  of  the  present  effort: 


i.  Using  the  transient  model,  the  effect  of  variations  in  the  cross-sectional 
shape  of  the  device  should  be  investigated.  The  present  shape  should  be 
analyzed  to  determine  how  accurately  it  matches  the  shape  as  specified  by 
the  manufacturer.  In  addition,  the  cross-sectional  shape  should  be 
optimized  with  respect  to  maximizing  the  heat  transport  capacity  and  the 
combined  effect  of  the  wetting  angle  and  the  shape  determined.  Shapes 
other  than  the  present  trapezoidal  cross-section  should  also  be 
investigated,  including  triangular  and  irregular  polygons. 

ii.  Again  using  the  transient  model,  the  effect  of  variations  in  the  amount  of 
working  fluid  should  be  studied.  The  results  of  this  investigation 
indicate  that  since  the  liquid-vapor  interface  changes  continually  along 
the  pipe,  and  hence,  micro  heat  pipes  are  very  sensitive  to  the  amount  and 
type  of  working  fluid  present.  In  addition,  the  wetting  angle  has  been 
identified  as  one  of  the  critical  parameters  in  the  determination  of  the 
maximum  heat  transport  capacity.  For  these  reasons,  determination  of  the 
effects  of  both  over  and  under  filling  on  the  transport  characteristics, 
along  with  the  development  of  a  method  by  which  the  amount  of  fluid  charge 
can  be  optimized  is  critical  to  the  successful  application  of  these 
devices. 

iii.  One  of  the  most  significant  applications  of  micro  heat  pipes  is  the 
incorporation  of  these  devices  into  silicon  or  gallium  arsenide 
semiconductor  chips  as  described  previously.  In  order  to  accomplish  this, 
information  is  required  to  determine  how  continued  reductions  in  the  size 
of  these  devices  affects  their  performance.  This  could  be  done  using  the 
previously  developed  analytical  model. 
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iv.  To  incorporate  micro  heat  pipe  technology  into  the  applications  mentioned 
previously,  we  should  understand  their  transient  behavior.  The  transient 
model  has  been  shown  to  accurately  predict  the  steady-state  behavior  of 
micro  heat  pipes,  but  the  predicted  transient  behavior  has  not  as  yet  been 
analyzed  or  compared  with  the  results  of  experimental  tests.  Therefore,  the 
experimental  investigation  should  be  continued  to  measure  the  transient 
characteristics  of  the  test  pipes.  Because  of  the  small  characteristic 
dimensions  of  the  test  articles,  the  thermal  measurements  would  need  to  be 
made  using  an  infrared  camera.  Upon  completion  of  the  experimental  tests, 
the  results  should  be  compared  with  the  transient  model  to  determine  how 
accurately  the  transient  behavior  effects  can  be  predicted.  Once  verified, 
the  transient  behavior  of  the  optimized  heat  pipe  could  be  evaluated. 
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VI.  NOMENCLATURE 


A  Cross-section  area 

a  Included  angle  of  the  groove 

B  Correction  factor  (constant) 

C  Constant 

Cp  Specific  heat 

d  Dimension,  diameter 

f  Friction 

g  Gravity 

j  Evaporation 

permeability 
K  Shape  factor 

L  Latent  heat,  length 

Mass  flowrate 
M  Molecular  weight 

P  Pressure 

Q  Heat  transfer  weight 

R  Universal  gas  constant 

r  Radius  of  curvature 

Re  Reynold's  number 

r^  Hydraulic  radius 

T  Temperature 

V  Velocity 

Vol  Volume 

w  Groove  width,  vapor  velocity 

wp  Wetted  perimeter 

W  Gas  velocity 

We  Weber  number 

Greek  Alphabet: 

a  Half  included  angle 

P  Dimensionless  variable 

S  Specific  heat  ratio 

e  Wick  porosity 

9  Contact  Angle 

U  Viscosity 

p  Density 

a  Liquid  surface  tension 

x  Shear  stress 

w  Groove  width 

Tilt  angle 

Subscripts; 

a  Adiabatic 

b  Boundary,  boiling 

c  Capillary 

d  Diagonal 

e  Evaporator,  entrainment 

ec  Evaporation-condensation 

eff  Effective 

f  Friction 


g  Groove 

h  Hydraulic 

i  Inner 

in  Input 

1  Liquid 

m  Momentum,  mean 
n  Nucleation 

o  Outer 

out  Output 

p  Pipe 

s  Solid 

v  Vapor,  viscous 
I  I  Parallel 
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